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ABSTRACT 

We describe our attempt to determine if gamma-ray burst (GRB) and afterglow emis- 
sions could both arise in external shocks for simple GRBs - bursts consisting of just 
a few peaks in their lightcurves. We calculate peak flux and peak frequency during 
the gamma-ray burst for ten well observed bursts using the same set of parameters 
that are determined from modeling afterglow emissions. We find the gamma-ray emis- 
sion properties for 970508 (which had a single peak lightcurve) fit nicely with the 
extrapolation of its afterglow data, and therefore this burst was likely produced in 
the external shock. One can explain two other bursts in this sample as forward shock 
synchrotron emission provided that the magnetic field parameter during the burst 
is close to equipartition, and larger by a factor ~ 10 2 than the afterglow value at 
~ lday. The remaining seven bursts cannot be explained in the external shock model 
even if we allow the energy fraction in electrons and magnetic field and the density of 
the surrounding medium to take on any physically permitted value; the peak of the 
spectrum is above the cooling frequency, therefore the peak flux is independent of the 
latter of these two parameters, and is too small by about an order of magnitude than 
the observed values. We have also considered inverse- Compton scattering in forward 
and reverse shock regions and find that it can explain the 7-ray emission for a few 
bursts, but requires the density to be 1-2 orders of magnitude larger than a typical 
Wolf-Rayet star wind and much larger than permitted by late afterglow observations. 

We have also calculated emission from the reverse shock for these ten bursts and 
find the flux in the optical band for more than half of these bursts to be between 
9th and 12th magnitude at the deceleration time if the reverse shock microphysics 
parameters are same as those found from afterglow modeling and the deceleration 
time is of order the burst duration. However, the cooling frequency in the reverse 
shock for most of these bursts is below the optical band, and therefore the observed 
flux decays rapidly with time (as ~ t~ 3 ) and is unobservable after a few deceleration 
times. It is also possible that the deceleration time is much larger than burst duration 
in which case we expect weak reverse shock emission. 

Key words: gamma-rays: bursts, theory, methods: analytical - radiation mecha- 
nisms: non-thermal - shock waves 



1 INTRODUCTION 

The localization of Gamma- Ray Bursts (GRBs) and the dis- 
covery of their X-ray afterglows by the BeppoSAX satellite 
in 1997 has greatly improved our understanding of GRBs 
over the last 7 years. These x-ray, optical, and radio after- 
glows are thought to be produced when an external shock 
heats the surrounding medium, with radiation being pro- 
duced via synchrotron from the heated material. We now 
also know from the spectroscopic confirmation of the GRB 
030329/SN2003dh connection (Matheson et al. 2003, Stanek 
et al. 2003) that at least some long duration GRBs are pro- 



duced by the collapse of massive stars. There is, however, 
considerable uncertainty surrounding the nature of the in- 
ner engine of GRBs, and we lack a definitive understanding 
for how the 7-ray emission is produced. This paper is an 
attempt to understand how 7-ray emission is produced in 
GRBs. 

Multiwavelength afterglow data have enabled us to do 
broadband modeling of late-time afterglows. This broad- 
band modeling results in the determination of burst energy, 
microphysical shock parameters, beaming angle, and envi- 
ronmental properties (surrounding density and stratifica- 
tion). Further improvements to our understanding of GRBs 
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requires analysis/modeling of both the GRB and afterglow 
together, which we undertake here. 

In this paper, we use parameters we determine for 10 
bursts by modeling their broadband afterglow emissions to 
extrapolate the radiation calculation back to the burst du- 
ration, with the goal of determining whether synchrotron 
emission from the forward shock can account for both the 
GRB prompt emission and the late-time afterglow. This is 
especially applicable for the ~10% of bursts with a sin- 
gle pulse FRED (fast rise, exponential decline) GRB light- 
curve, where a single external shock is expected to produce 
the emission. This exercise is, however, carried out for all 
bursts in our sample, including those with moderately com- 
plex GRB lightcurves. 

Internal shocks were suggested as a mechanism for pro- 
ducing 7-ray emission because external shocks are not ca- 
pable of producing rapid variability seen in many GRB 
lightcurves, whereas variability arises naturally in internal 
shock models, reflecting fluctuations associated with the 
central engine (Rees & Meszaros 1994; Piran, Shemi & 
Narayan 1993; Katz 1994). For GRB lightcurves consisting 
of a single peak or just a few peaks, this rationale for internal 
shocks does not apply and such bursts could be produced in 
external shocks. 

The determination of kinetic energy release in relativis- 
tic ejecta for ten bursts by modeling their broadband after- 
glow lightcurves suggests that the efficiency for 7-ray pro- 
duction is typically in excess of 50% (Panaitescu & Ku- 
mar, 2002 hereafter PK02). Such a high efficiency cannot 
be achieved in internal shocks; some published claims to the 
contrary (e.g. Beloborodov, 2000) achieved high efficiency 
by colliding shells with very large relative Lorentz factor 
(hereafter LF), however in this case the emergent spectrum 
peaks at energies much larger than observed values. Exter- 
nal shocks, on the other hand, can very efficiently convert 
bulk kinetic energy to radiation. 

In addition to the problem of efficiency for the internal 
shock model we describe below other reasons for considering 
the external shock model for the generation of 7-ray emission 
for many of the ten bursts we consider in this paper (table 

1 lists the ten bursts). 

The 320-1090 kev light-curve for 990123 consisted of 
two broad peaks of duration ~ 10s each. Comparing this 
time scale with the deceleration time (td) of <50s - which 
is inferred from the peak of the prompt optical emission 
- suggests that 7-ray emission is produced within a factor 

2 of the deceleration radius 1 . In the internal shock model 
for 7-ray production, the near equality of the radius where 
shells collide to produce 7-rays and the deceleration radius 
is a coincidence, whereas in the external shock model this 
is what one expects. It should be noted that the short time 
scale variability seen in 990123 (Fenimore et al. 1999) had 
an amplitude of ~20% and could have arisen due to small 
scale turbulence in the shocked fluid. The observed low en- 
ergy spectral index a for this burst was 0.4 (/„ oc v aA ) 



When the outermost 7-ray producing shell undergoes decelera- 
tion and is heated by the reverse shock it produces optical flash, 
and its radius increases as ~ t 1 / 4 for t><cj/4 (where t is the ob- 
server time). Therefore, the increase in the radius for a 5- fold 
increase in time is less than a factor 2. 



whereas in internal shock models we expect a ~ —0.5 due 
to short cooling time for electrons or low cooling frequency 
(Ghisellini et al. 1999). 

The lightcurve for GRB 970508 was a FRED, 980519 
was similar to a FRED, and 000301c lightcurve was per- 
haps a FRED (Smith et al. 2002), however because of the 
low temporal resolution of the Ulysses observation (0.5 sec) 
we are unsure of it. One might expect these bursts to arise 
in an external shock. Two other bursts in our sample of 
ten - 980703 &: 991208 - had lightcurves consisting of two 
smooth peaks, and therefore are good candidates for a pos- 
sible origin in an external shock. GRBs 990510 & 991216 
lightcurves had more fluctuations than the bursts mentioned 
above, however they each had two broad peaks and a num- 
ber of sub-pulses superimposed on them, and do not require 
internal shock to produce this modest variability. There are 
no lightcurves available for the remaining two bursts in our 
sample, 000418 & 000926, which were both detected by the 
IPN. It turns out that for all of these bursts, except 970508, 
the simplest theoretical model of synchrotron emission in 
the forward shock fails badly to explain their 7-ray emission 
(§2). Moreover, none of the possibilities we explore in the 
general framework of an external shock model seem to work 
satisfactorily. 

In two GRBs (990123 and 021211), a bright, steeply 
falling off (~ t~ 2 ) early optical emission was observed. This 
has been explained by radiation from the reverse shock 
heated ejecta from the explosion. We have seen this emission 
from only these two bursts, while there are many cases for 
upper limits within a few hundred seconds after the GRB 
time and even a few bursts (e.g. 030418 and 021004) with 
early afterglow detections that do not exhibit the bright, 
steep optical decay. In this paper, we also estimate the re- 
verse shock emission at deceleration for these ten bursts, and 
discuss possible reasons for numerous non-detections. 

§2 outlines the afterglow fitting and describes our 
method for calculating the flux and the peak frequency dur- 
ing the GRB. A comparison between the theoretical calcu- 
lation and 7-ray observations is also described in §2. Some 
alternate possibilities to explain the 7-ray observations such 
as inverse Compton in the forward or the reverse shocks, 
pair enriched ejecta, or high density clumps in the circum- 
stellar medium, are discussed in §3. In §4, we discuss reverse 
shock optical emission at deceleration. 



2 AFTERGLOW TO 7-RAY EMISSION 

The afterglow modeling is described in detail in Panaitescu 
& Kumar (2001) and (2002). Briefly, we determine the colli- 
mated fireball dynamics by numerical integration of a simpli- 
fied set of jet propagation equations, keeping track of radia- 
tive loss of energy due to synchrotron and inverse-Compton 
emissions. The synchrotron peak and cooling frequencies are 
calculated by assuming that a certain constant fraction of 
the thermal energy of the shocked fluid is imparted to elec- 
trons and magnetic field. The effect of IC loss including the 
proper Klein-Nishina cross-section is included in the calcu- 
lation of the cooling frequency. The observed lightcurves are 
calculated by integrating the emissivity over equal arrival 
time surface. All of the unknown parameters, which include 
jet opening angle, the total energy release in the explosion 
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(which is is the sum of the kinetic energy given in PK02 and 
the energy in 7-ray radiation), the fraction of energy in elec- 
trons (e e /), and the fraction in magnetic field (es/), are ob- 
tained by fitting the observed light-curves and the spectrum 
with the theoretically calculated curves by a % 2 minimiza- 
tion. The parameter e' B f, which determines the minimum 
thermal Lorentz factor of electrons, is 0.1 for all bursts for 
which p > 2 (PK02). Since the high energy spectral index 
during the burst gives p > 2 we set e' e f =0.1 for all bursts 
in our calculations during the gamma-ray burst. 

Using these parameters we estimate the frequency 
where the spectrum (vf„) peaks and the flux (/„) at this 
peak at deceleration (which we assume is half of the 7-ray 
burst duration). The results for ten bursts are summarized 
in Table 1 for a uniform circumburst density (wind circum- 
burst medium had similar results, and for brevity are not 
listed here). The theoretical results are compared with the 
observed data for these bursts (see Table 1). Note that for six 
out of ten bursts in the table the peak frequency during the 
burst is within a factor of about 2 of the the observed value 
which we consider a reasonably good agreement. However, 
in four of these cases the theoretical peak flux is smaller than 
the observed value by an order of magnitude or more. For 
GRB 970508, which was a single peaked FRED burst, the 
fluxes are in good agreement. Therefore, for this burst the 
7-ray emission could arise in an external shock; it is highly 
encouraging to see the forward shock model works so well 
to explain observations all the way from the 7-ray emission 
at 10s to radio at 100s of days. However, for 000301C which 
was also likely a single pulse FRED, and 980703 & 991208 
each of which contain two simple peaks in their 7-ray light- 
curve and are therefore good candidates for external shock 
mechanism for 7-ray production, the discrepancy between 
theory and observation is large. 

To understand how sensitive the 7-ray emission is to 
errors in afterglow modeling and parameter determination, 
and to consider some possible solutions within the frame- 
work of the external shock model, we present an analytical 
derivation of the main results. 

The forward shock synchrotron injection frequency, Vi / , 
and the flux at the peak of the F v spectrum are (Wijers & 
Galama, 1999) 



v if (t) = 
F P f(t) = 



Q.98qBj 2 T 

2irm e c{l + z) 

N e P Vp Y{l + z) 
4ndl 



(1) 
(2) 



where q & m e are electron charge and mass, m p is proton 
mass, cLl = 2c^/l + z[(l + z) 1 ^ 2 — l]/Ho is the luminosity 
distance, 74 = e' e f(m p /m e )(T — 1) is the minimum thermal 
LF of electrons (the electron distribution for 7 > ji is as- 
sumed to be a power-law of index p, i.e.dNe/d'y oc 7 _p ), 
N e = 'iirAR 3 ~ s mp 1 /(3 — s) is the total number of swept-up 
ISM electrons, po = AR~ a is the density of the medium just 
ahead of the shock, V is the bulk LF of shocked gas, 



B = 4cF [2Tre Bf AR~ s ] 



1/2 



t is the observer time, and 
_ 1.04q 3 B 



R= (4-s) C rV(l + 2), (3) 



(4) 



is the power radiated per electron per unit frequency, in the 
shell comoving frame, at the peak of the synchrotron spec- 
trum. The numerical factors of 1.04 in the above equation 
and 0.98 in equation Q are taken from Wijers and Galama 
(1999) for p = 2. 

The synchrotron injection frequency and peak flux, at 
deceleration, for the particular cases of s = & 2 are given 
below 
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F P f(t) = 
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s = 
s = 2 



(6) 



where A, — A/(5 x 10 11 ) g cm -1 , £ is the isotropic equiv- 
alent of energy release in the explosion, td is the observer 
frame deceleration time in seconds, and an integer subscript 
n on a variable X, X n , means X/10 n . In the derivation 
of the above equations we substituted for V using the equa- 
tion 4tt AR 3 ~ a c 2 (r 2 -l)/ '(3 -s) =£/2 at deceleration which 
states that half of the original kinetic energy of the explosion 
{£/2) is deposited into swept-up ISM; the LF at deceleration 
is given by: 



4.17 (^A 

\ no J 



1/8 



0.74 



(l + z)£52 



A,t r , 



(i+jO 



1/4 



3/8 



s = 
s = 2 



(7) 



We next calculate the electron cooling frequency. For 
this we need the Compton Y parameter, defined as Y = 
tt J d r y e y 2 (dn e /d'y e ), with tt being the column density 
of electrons times the Thomson cross-section. The Y- 
parameter is obtained by solving the equation describing 
the radiative loss of energy - 



d_ 



2 dn e 
m e c 7 e - — 

dje 



dn e (l + Y)a T B 2 j 2 c 
dj e 6tt 



(8) 



where t' is the comoving time, B is the magnetic field which 
we assume is uniform, and gt is the Thomson scattering 
cross-section. Consider the comoving frame down-stream 
fluid velocity to be v. This relates t' and the comoving ra- 
dial coordinate r' viz. dr' = vdt' . Changing the independent 
variable from t' to r' and integrating the above equation over 
the electron distribution we find 



d(e e U) 
dr' 



(1 



Y> T B' 2 cn e7 f 



6lTV 



(9) 



where U is the thermal energy density of shocked fluid, and 
7I is the average y 2 . Integration of this equation over r' for 



a highly relativistic shock (v 
(i>i v c ) we find 



Y(l + Y) 



4s, 



c) and highly radiative fluid 



(10) 



The calculation of the cooling LF of electrons, j c , at 
deceleration, is straightforward and is given below; 7 C is the 
LF of electrons that lose their energy in a time available 
since crossing the shock front averaged over the population, 
given by ^~t d /3 (2 - s)/2 . 
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7c (id) = 



37rm e c(l + z) 
2a T B'H c V d {l + Y) = 
3(l + z) 1 - s m e r^- 3 
64a T s B Ac(4ct d ) s t c (l + Y) 



(11) 



This in turn is used to calculate the synchrotron frequency, 
in observer frame, corresponding to the LF 7 C , and is referred 
to as the cooling frequency (v c ). 



Vc{t d ) 



2nm e c(l + z) 



° r = { 2.2 x IP' 7 J £ ^ )1/2 , 



(12) 



where we have made use of equations Q, 110L and 11H . We 
see that the cooling frequency at deceleration is typically 
much smaller than the synchrotron injection frequency (see 
cq. |5J and thus the peak of the vf v spectrum will generally 
be at Vif. For the case where Vif > i/ c , the flux at the peak 
of the vf v spectrum is obtained by using equations JSJ, Jfjjl 
& (I I 12 1 > and is given here 



F Vv {U) 



(g 52 t d ) i/2 (i+^r i/2 
' < f 4 2 i(i+^) i/2 -i] 2 

6 x 1CT 4 mjy s = 
3.5 x 10" 4 mjy s = 2 



(13) 



Note that the flux is independent of eb and the density of 
the ISM - the two parameters that have the largest error 
associated with them in afterglow modeling (PK02). Using 
equations J5J and 11311 we calculate the peak frequency and 
the flux during the 7-ray burst and the results are in good 
agreement with the numerical calculation result presented 
in table 1. 

Table 1 shows that for all those bursts for which the ob- 
served and theoretical peak frequencies agree within a fac- 
tor of two 2 (6 cases altogether), the theoretically calculated 
fluxes, with the exception of 970508, are too small by an 
order of magnitude compared to the observed fluxes. 

We see from equation 1131 that the flux at the peak of 
vf v depends only on £, which is a very well determined pa- 
rameter, and e' e f which is 0.1 for all the bursts with p > 2, 
and therefore there is no way to reconcile the difference be- 
tween the observed flux and the theoretical expectation in 
the simplest version of the synchrotron emission in the ex- 
ternal shock model. Any error in the parameter determina- 
tion from afterglow modeling does not affect our calculation 
of the 7-ray flux. In other words, even if we consider ebs 
and the ISM density n during the burst to take on values 
completely unrelated to what is determined from afterglow 
modeling the observed flux cannot be reconciled with the 
theoretical expectations in the forward shock model. This 
is a very robust result and effectively rules out synchrotron 
origin for 7-ray emission in the external shock model for 
these bursts. 

2 We have assumed that the peak frequency for 000301c was ~ 
500 keV, at the higher end of the peak frequency distribution, 
since it was a fairly short burst (8.4 seconds). 



In four cases, viz. 980519, 980703, 990123 & 000418, the 
theoretically calculated peak frequencies are much smaller 
than the observed value; these are the four cases with the 
smallest ebs as determined by the afterglow modeling (see 
Table 1). Could a larger s B f at early times, as in the case 
of 021211 (Kumar & Panaitescu, 2003, hereafter KP03), ex- 
plain the peak frequencies for these four cases? If £s/ were 
to be 0.5 during the burst for 980519 & 000418 we can ex- 
plain the 7-ray emission for these GRBs. However, even 
if we set e B f = 1 during the burst for 980703 & 990123 
the synchrotron frequency falls short of the observed value. 
Could the gamma-ray burst in these cases arise as a result 
of inverse-Compton (IC) scattering of the synchrotron radi- 
ation in the forward or the reverse shock? We consider this 
possibility, and some others, in the next section. We also in- 
vestigate whether the peak flux of the IC component might 
be able to match the observed flux for the other five bursts 
which have too small synchrotron flux. 



3 GAMMA-RAYS IN EXTERNAL SHOCK: 
SOME ALTERNATE POSSIBILITIES 

We consider below (§3.1) a combination of synchrotron and 
inverse-Compton processes in the forward and reverse shocks 
to determine if this could explain the 7-ray emission prop- 
erties for the nine "problem bursts" in our sample of ten 
discussed in the last section. In §3.2 we discuss if a collision 
between the GRB ejecta and a high density clump might 
be able to explain the large 7-ray flux for the five of the 
bursts in our sample, and in §3.3 we look into the effect of 
electron-positron pair loaded ejecta on 7-ray emission. 

3.1 Inverse Compton in the external shock 

We investigate the effect of IC in external shocks - forward 
as well as the reverse shock - to see if the observed 7-ray 
emission for the bursts in our sample could be explained by 
the IC process. 

Consider the flux at the peak of the synchrotron ra- 
diation vf v spectrum, v p , to be f v . We consider inverse- 
Compton scattering by a population of electrons that could 
be distinct from the population that gives rise to the syn- 
chrotron radiation. For instance, the synchrotron emission 
could be produced in the reverse shock (RS) and the IC 
scattering in the FS. Let us take the minimum thermal 
LF of electrons in the IC-scattering region to be 7 m i n = 
min(7i,7o), the electron distribution to have a break at 
76 = max(7i,7 c ) such that for 7 > 75 the electron distri- 



bution is proportional to 7 



-P-i. 



7» 



7 C are the injection 



and cooling LFs for electrons. The peak of the IC radiation 
(for vf v ) is at 



ic 



(14) 



If the optical depth of the medium to Thomson scattering 
is tt, then the flux at the peak for the case where 74 <C j c 
is 



fu p ~ T T (j m inhbY fu p 



(15) 



For 7 C < 7,:, there is a slightly different relationship. 

The optical depth at deceleration in the forward shock 
is given by 
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5 



Burst 


£ 

(xlO 52 ergs) 


£ Bf 




fj)(td) 
(keV) 


, obs 

(keV) 


iv p (td) 
(mjy) 


J Up ,obs 

(mJy) 


Refs. 


970508 


7.92 


4.5 x 10" 


2 


98.4 


100 


0.51 


0.6 


2,4 


980519° 


99.8 


3.5 x 10- 


-5 


7.0 


700 


1.25 


0.3 


3,4 


980703 


0.40 


3.0 x 10" 


3 


0.8 


370 


0.21 


0.3 


3 


990123 


277.5 


9.9 x 10- 


-4 


6.6 


780 


1.24 


9.0 


1,2,3,5 


990510 


55.5 


5.0 x 10- 


3 


184.9 


160 


0.27 


4.0 


1,2,3 


991208 


11.9 


3.8 x 10- 


2 


118.7 


190 


0.78 


20.0 


1,3,6 


991216 


100.0 


6.2 x 10- 


3 


153.4 


410 


1.0 


30.0 


1,3 


000301c 


8.19 


1.0 x 10- 


1 


541.6 


500 fe 


0.04 


0.3 


1,3,7 


000418 


8.05 


1.6 x 10- 


2 


25.2 


280 


0.32 


0.1 


1,3 


000926 


39.7 


8.1 x 10- 


2 


198.4 


130 


0.21 


0.7 


1 



0.1, e. 



0.5 for all bursts, no, 9q are equal to values found from afterglow modeling. 



References: (1) Mazets (personal communication); (2) Amati, L. et al. 2002; (3) Jiminez, R. et 
al. 2001; (4) Nicastro, L. et al 1999; (5) Briggs, M. S. et al 1999; (6) Hurley, K. et al 2000; (7) 
Smith, D. A. et al 2002 



a Redshift not known for this burst, 



1 used. 



No v. 



p,obs 



available for this burst; a value of 500 kev assumed. 
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<Tt£ 



4nm p c 2 ^ 2 l R 2 d , 
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where Rd = (4 — s)ctdT d /(l + z) is the deceleration radius, 
Yd is the LF at deceleration, and Fo ~ 1.51^ is the initial 
LF of the ejecta (see eq. 1191 . Using equation J7J this can be 
rewritten as follows 



TT 



3.9 x l(r 9 nj| /4 £^ 4 [td/(l + z)] 1/4 s = 
8.7 x 10- 4 ^ /2 f 5 - 2 1/2 [(l + z)/t d } 1/2 s = 2 



(17) 



We see that the optical depth for a uniform density ISM is 
very small, and therefore the inverse Compton flux due to 
scattering in the forward shock region, for s — 0, is likely to 
be too small to be observationally interesting. 

The optical depth to Thomson scattering of the ejecta 
at deceleration can be obtained directly from equation 11711 
by recognizing that the mass of the ejecta is larger than the 
swept-up ISM mass by a factor F d . Thus, 



TT = 



1.5 x 10- 6 n B /B £U S itd/(l + zT l/s 



5.4 x 10- 2 ^ /4 £- 1/4 [(l + z)/t d ] 



3/4 



(18) 



3.1.1 Reverse Shock Break Frequencies and Peak Flux 

To complete the calculation of inverse Compton scattering 
of synchrotron emission produced in the reverse shock (RS) 
region we provide below the synchrotron characteristic fre- 
quency and the flux from the RS (see KP03 for details). 

The thermal energy per proton in the RS at decelera- 
tion, e p , can be calculated using the following pair of equa- 
tions 

9 -| -1/2 

I - 2 ( — - - j . (19) 



2 Vr + r d 



El 

r„ 



which fits the results of numerical calculations to better than 
8% in the Newtonian, Relativistic and intermediate regimes; 
where n ej is the comoving density of the ejecta, no is the 
ISM density, and To is the initial LF of the ejecta. 



oej/no at the time when the re- 



It can be shown that n e 
verse shock arrives at the back end of the ejecta (which is 
approximately equal to the deceleration time for the ejecta) 
is 1.5F 2 , for a uniform density ISM and 3.5F 2 , for s — 2 
medium. The reverse shock in this case is neither highly rel- 
ativistic nor Newtonian. Using the above equation we find 
the thermal energy per proton in the RS in this case to be 
0.13m p c 2 (0.067m p c 2 ) for a s = (s = 2) medium. 

The injection frequency at deceleration for RS is smaller 
than the FS by a factor of r d /0.13 2 for uniform ISM and is 
given below 



Vir (td) — 



0. 07 qm p 



5/2^1/2^/2 „l/2^- S /2 p2 



-Br ' 



'•() 



(2tt) 1/2 ml(l + z) 



(20) 



This equation, and a similar one for s = 2, can be rewritten 

as 

f,\_ 1/2 12 
fir \td) — ^ Br 

37nJ /4 f 1 2 /4 (l + ^)- 1/4 t; 3/4 ev s = 



610 Al^t' 1 ev 



s = 2 



(21) 



where we have used equation @ to eliminate Td - The cooling 
frequency in the reverse shock region, when vw > v C r, is 
given by equation 111 2ti with appropriate values of eb and 
£g for the reverse shock. However, the reverse shock v cr is 
typically larger than Vi T (see Table 2), so equation HUH is 
not a valid approximation for the Compton Y parameter 
any longer and we must also use an appropriately modified 
form version of equation 1121 . 

The flux at the peak of the reverse shock /„ spectrum 
at deceleration is larger than the peak flux from the FS by 
a factor Yd and can be written as 



Fpr (td) 



(SSBrA^q'Sjl + zy/ 2 

1 (Actd) s/2 [(3 - a) : 



m e mpC 3 d'£T d 



,1/2 



(22) 
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Fpr (td) 



.1/2 



[(l + *)i/2_i]a 

J 6.8 x 

X \ 2.7 x 10 5 (1 + z)'^ 4 £U 4 A 3 J 4 t- 3/4 mJy s = 2 



xl0 2 (l + 2 ) 3 /% 9 / 8 nf/V /8 m Jy s = (23) 



The self-absorption frequency in the reverse shock region, 
VAr, is often as large as the cooling and the injection fre- 
quencies, and therefore should be taken into consideration 
in the calculation of observed flux. The self-absorption and 
cooling frequencies and the Compton Y parameter need to 
be calculated together in a self-consistent way (as we do for 
all of our numerical calculations). However, when electron 
cooling is dominated by the inverse-Compton scattering and 
max{i/ lr ,i/ cr } > VAr, the calculation of self-absorption fre- 
quency is considerably simplified and is given by (KP03) 



8.4 x lO- 2 ^ W /8 
49^/Vev 



(1 + *) 



-3/8 



ev s = 
s = 2 



(24) 



where a depends on the ordering of v ir , v c and VAr, and is 
equal to 1/3 if v c > v ir > vat, and — p/2 when v c > vat > 



3.1.2 Inverse Compton results 

The peak frequency and flux for the IC radiation is cal- 
culated as described in §3.1. The calculation of the syn- 
chrotron injection frequency is straightforward & is carried 
out as described in §2 for the FS and §3.1.1 for the RS. 
The synchrotron self-absorption frequency is typically small 
in the FS and is unimportant for IC calculation. However, 
in the RS the self-absorption can be larger than the cool- 
ing frequency and these frequencies must be calculated self- 
consistently; we calculate these frequencies numerically. 

There are four cases of the IC scattering to consider: 
synchrotron in the FS and IC in either the FS or the RS, 
synchrotron in RS and IC in the RS or the FS. We have 
investigated these cases numerically, and we have explored 
the parameter space - £, EBr, e'en £ Bf, s' e f, n - for each burst 
to determine if the observed 7-ray peak frequency and flux 
could be explained by the IC radiation, either for a uniform 
density circum-burst medium or a wind-like medium (p oc 
R~ 2 ). The results for each burst are described below. 

For a uniform density ISM the synchrotron-IC mecha- 
nism in external shock offers a vanishing parameter space 
that is consistent with the gamma-ray emission properties 
for GRB 990123. However, for a s — 2 medium we find some 
solutions where the synchrotron emission produced in the 
forward shock undergoes inverse-Compton scattering in the 
reverse-shock region. The density required for these solu- 
tions is ~ 10 2 times that normally associated with Wolf- 
Rayet star winds, and much greater than what is found 
from modeling of early and late time afterglow observa- 
tions. Other parameters such as the energy in the explosion, 
micro-physics shock parameters in the forward and the re- 
verse shock are roughly consistent with the afterglow ob- 
servations. However, the low energy spectral index for the 
7-ray spectrum (a) is -0.5 whereas the observed index is 0.4. 
Therefore, we do not have a fully self consistent solution for 



the gamma-ray emission properties for 990123 in the exter- 
nal shock model. This is surprising in the light of the argu- 
ments for external shock (see §1) for this burst, and perhaps 
suggests that there may be some other mechanism produc- 
ing the 7-ray photons that is completely different from the 
standard internal/external shocks model. In the next sub- 
section we explore if high density gas near the deceleration 
radius could explain the 7-ray emission. 

In the case of 980519 there are IC solutions for s = 2 
medium where the synchrotron radiation is produced in the 
forward shock & IC in the reverse-shock, and the parameter 
space consists of A, in the range of 15 to 100 which is at least 
a factor of a few larger than the value of A, ~ 3.5 determined 
from afterglow modeling. However, the IC solution requires 
eb ~ 10~ 4 in the forward shock that is much smaller than 
the value of 0.1 we find for this burst from the afterglow 
fitting when s = 2. The RS optical peak flux of ~ 11th 
magnitude is perhaps not a problem for this solution. 

For 980703 there are solutions found for s = 0. How- 
ever, these solutions require £ > 10 55 erg and the RS op- 
tical flux is larger than 1 Jy, and therefore these are not 
acceptable solutions. The solutions we find for this burst 
with s = 2, which involve synchrotron in the FS and IC in 
the RS, require A„>5 and other parameters are roughly con- 
sistent with the values we find from the afterglow modeling; 
the optical flux from the RS is ~ 20 mJy. If we ignore the 
somewhat high density requirement, this burst could per- 
haps be produced as IC in the external shock. 

For 990510 no solution is found that is in agreement 
with the observed properties of this burst. The same is 
true for 991208, 991216. To be precise, there are no so- 
lutions found when the density of the medium is taken 
to be uniform. However, for a pre-ejected wind medium 
there are regions in the multidimensional parameter space 
(£,A,,eB,s' e ) that give gamma-ray flux and peak frequency 
in agreement with observations for these bursts where syn- 
chrotron emission is produced in the forward shock and the 
inverse-Compton scattering takes place in the reverse shock 
region. The problem is that for all of these "solutions" A, 
is greater than about 10 2 which is larger than what we ob- 
tain from afterglow modeling by two orders of magnitude, 
and too large for winds from Wolf-Rayet stars. Moreover, 
the large ebt in the reverse shock for these solutions gives 
rise to optical R-band flux of about 10 Jy, or 6th magnitude, 
which is unlikely to have gone unnoticed. Therefore we do 
not consider these solutions physically acceptable. 

We find two different IC "solutions" for 000418 for a s = 
2 circum-burst medium; one of which is synchrotron in the 
FS and IC in the RS, and the other is synchrotron in the RS 
and IC in the FS. The first scenario requires A*>25 (which 
is larger by factor ~ 10 2 than determined from afterglow 
observations) and £s/ is between ~ 10 -5 and 2x 10~ 4 which 
is smaller by two orders of magnitude than the value we find 
from afterglow modeling. The peak optical flux from the RS 
is between 20 mJy and 3 Jy; the lower flux value in this case 
certainly poses no difficulty with observations. The latter 
scenario requires A»>400 and es r ~ 10~ 5 which are unlikely 
to be realized in nature. 

There are also two solutions for 000301c with s = 2. For 
FS synchrotron & IC scattering in the RS, A* is about an 
order of magnitude larger than that found by afterglow mod- 
eling, but the optical flux from RS is not a problem, being 
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<0.5 Jy. The magnetic field parameter in the forward shock 
required for the IC solution is, however, several orders of 
magnitude smaller compared with the value we find from af- 
terglow modeling in a wind-like medium for this burst. This 
together with the required high density for the CBM makes 
this solution unacceptable. For the case of RS synchrotron 
& IC scattering in the FS, the allowed parameter space to 
explain the 7-ray observations require the energy to be a 
factor of 10 smaller than the observed value, and A*>100. 
The magnetic field parameter is also 100 times smaller than 
the afterglow value, so we do not consider these solutions 
viable. For a uniform density medium the IC solution for 
000301c requires £ > 10 55 which is more than two orders 
of magnitude larger than the energy determined from either 
the 7-ray fluence or the afterglow emission. The situation 
for 000926 is a bit worse than 000301c. 

We also note that the IC flux for 970508, in any of the 
four possible combinations, at the observed peak is too small 
to be a significant contributor to the observed flux. 



3.2 Effect of density clumps in the ISM on 
gamma-ray flux 

In this subsection we investigate whether a dense clump of 
gas in the circum-burst medium (CBM) might increase the 
flux in the 7-ray band and thereby explain flux observations 
at the peak of the vf v spectrum for some of the bursts in 
our sample. 

Consider a dense clump of angular size greater than 
r^ 1 , and proton number density n; Fo is the initial Lorentz 
factor (LF) of the ejecta which we assume does not decrease 
until it runs into the clump. For the calculation of radiation 
such a clump can be treated as a spherical object. We take 
the external density to be sufficiently high that the forward 
shock LF is less than To and the reverse shock is relativistic; 
the parameter £ = Y^n/n e j > 1 determines the thermal LF 
of protons in the forward and the reverse shocks; n e j is the 
density of the ejecta when it hits the clump. For large n the 
thermal LF of protons in the FS is 7 Pi/ = r £~ 1/4 /2 1/2 and 
in the RS it is £ 1/4 /2 1/2 . 

The thermal energy density in these shock regions 
is 4nm p c 2 7p ^ ~ nT 2 m p c 2 ^ 1 '" 2 , and therefore the mag- 
netic field strength B oc (e s/ n) 1/2 r £~ 1/4 . The syn- 
chrotron injection frequency in the forward shock is Vi oc 

The number of swept up protons at deceleration is 
obtained by equating the thermal energy of protons (in 
lab frame) with half the energy in the explosion i.e., 
JV p m p c 2 7 2 i/ = £/2 or N p = ££ 1/2 /(m p c 2 r 2 ). 

Let us assume that the distance of the clump from the 
center of the explosion is R c and the forward shock trav- 
els a distance of SR C before the shocked material acquires 
half the energy of the explosion. Therefore, N p = AimR^S = 
0). The GRB duration (in the observer frame), 
if the 7-rays are produced due the ejecta colliding with 
the dense clump, is t 7 ~ Rc/{ c l P ,f) ~ RcS, 1 ^ 2 / {cYq). 
Combining this with the equation for N p we find: tlTq ~ 
£ 2 £ /(47rm p c 5 i^<5). Substituting this back into the equation 
for injection frequency we find Vi oc s 1 J^£ , e f 2 £ ^i" 3 ^ -1 / 2 . 

The time scale in the lab frame for electrons to cool 
in the forward shock is ~ SR c /(cy p j), from which we cal- 



culate the cooling frequency to be v c oc (nes/) 3 ^ 2 (1 + 
Y)- 2 (8R C )- 2 oc £ B f 2 (t 7 £ )- 1/2 rr 1 5- 3 / 2 (l + Y)- 2 . 

The flux at the peak of the synchrotron spectrum is 
f p oc £ 1 JJn 1 ^ 2 £ . For v c < i>ij, expected for a high den- 
sity clump, the Compton Y-parameter is ~ (e e f /EBf) 1 '" 2 , 
and the flux at the peak of vf v i.e. at Vif is proportional 
to £ 1 ^ 2 tl/ 2 S~ 1 ^ 2 e~^ ^ 2 s' e f 1 • Note a weak dependence of the 
peak flux on S = SR C /R C . 

The distance 5R C the FS moves before the GRB ejecta 
is decelerated is obtained by calculating the density of the 
ejecta at R c \ n e j ~ f /[47ri? 2 m p c 2 Fo max(i»c, R c /Fq)], where 
t* is the duration of the central engine in the lab frame. 
Using this expression and the definition of £ we find that 
S ~ max(£*/i 7 ,£~ 1/2 ). 

Since the peak flux is proportional to (Y 1 / 2 , a factor of 
ten increase in the flux requires £>10 4 or the clump density 
(lO^F 2 ,). Therefore, for r ~ 10 2 , n is of order 7i e j 
and we find that in order to explain the gamma-ray flux in 
the clump-ejecta collision the density of the clump needs to 
be similar to the ejecta as in the internal shock model! Very 
bright early afterglow will be produced in such a collision 
which might pose a problem for this scenario. We note that 
the increase in the flux is accompanied by an increase in the 
peak frequency, both of which are proportional to <5~ 1//2 , but 
the latter quantity can be easily adjusted by a decrease in 
EBf or e' e to match the observations. 

3.3 Effect of e pairs present in the ejecta 

In this subsection we consider whether e pairs present in 
the ejecta can make inverse-Compton radiation match the 
peak flux and frequency at deceleration. Pairs would soften 
the reverse shock spectrum, and the larger optical depth 
could increase IC flux at the peak. Adding N± pairs per 
proton to the ejecta increases the number of radiating par- 
ticles in the reverse shock region thus lowering the energy 
per particle. This decreases Vi r by a factor of N± while in- 
creasing the peak flux of the f v spectrum by a factor of 
N±. The cooling frequency for the highly radiative regime 
is not affected, nor is it changed when the Compton-Y pa- 
rameter is much less than one; it is, however, affected when 
Vi r < v cr (more common in the reverse shock emission than 
Vi r > ^cr, see Table 2) and Y > 1, increasing it by a factor 

jy2(p-2)/(4-p)_ 

The flux at the cooling frequency, which is the peak of 
the reverse shock synchrotron vf v spectrum, is proportional 
to N^ 2 ~ p M( 4 p ' . We see that for p — 2 the cooling frequency 
and the peak flux are independent of N±. The peak flux 
increases by a factor of TV 3 / 5 for p = 1.5 (decreasing for p > 
2). For inverse-Compton scattering of these reverse shock 
synchrotron photons by the ejecta, the IC peak, v p> ic = 
j 2 v c , will increase by a factor of 7\r^ p ~ 2 ' ) ^ 4 ~ p ' 1 . We see again 
that Vpjc is not affected when p — 2, but can go up (p > 2) 
or down (p < 2) for other values. The IC flux at the peak 
changes by a factor of jV^ 2 ~ p ' )/ ' ( - 4 ~ p ' ) (from its value without 
pairs present), which, similarly, does not change much for 
values of p around 2. The IC peak flux decreases for p > 2 
and at most, can be increased by a factor of ~ N± when the 
electron distribution is very hard, i.e. p<1.5. 

It can be shown that the peak flux for synchrotron pro- 
duced in the FS and IC in the slow cooling RS (with Y >• 1) 
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is proportional to jV^( 2_p )/( 4_p ) ; decreasing for p > 2. If the 
RS was in the highly radiative regime, the IC flux would 
increase as N±; however, the injection frequency decreases 
rapidly as N± , whereas v cr is independent of N± in the fast 
cooling regime. We therefore expect Ui r to become less than 
v cr as N± becomes larger than a certain value and the slow 
cooling regime considered earlier once again applies. The re- 
verse process - synchrotron in the RS and IC in the FS - 
has also lower peak flux for p > 2. Thus, we find that pairs 
present in the ejecta are not likely to be able to account for 
the theoretical IC and observed 7-ray flux difference. 



4 EARLY AFTERGLOW EMISSION 

It is generally believed that the steeply falling off early after- 
glow emission observed from GRBs 990123 and 021211 was 
produced by the reverse shock heated ejecta from the explo- 
sion. This emission falls off roughly as i -1 ' 7 and flux falls 
below the forward shock emission level after about 10 — 20 
minutes. We use the equations in §3.3 to calculate the ob- 
served flux in the optical R-band at deceleration for this 
sample of 10 bursts. For the case of v ir < v ar < vr, the flux 
at v r ~ 2 eV, the R-band in observer frame, is given by 



fit (td) 



[(1 + 2)V2 _ 1]2 



400 x 3.3- 



d,l 



mjy 



I 584 x 30 — £52^, 4 t d l mjy 
Whereas for the case where Vi r < vr < v cr the flux is 

fn (td) 



s = 
s = 2 



U/2 



p-1 4-p £±8 2+1 -IE 

290 x 3.3— (1 + z) — £ 52 s n 8 t dl s mjy 

p-1 3 P + 2 2p+l 

4.8xl0 4 x 30— £ r %A, 4 t. 



(l + z) 4 mjy 



Table 2. Predicted reverse shock flux using afterglow param- 
eters for homogeneous external medium 



Burst 


R a 


(eV) 




v cr 
(eV) 


(eV) 




970508 


9.8 


1.0x10" 


-2 


5.8 


3.9 x 10' 


-2 


980519 6 


17.1 


3.8 x 10 


-4 


5.1 x 10 5 


2.8 x 10" 


-3 


980703 


17.4 


2.4 x 10 


-5 


6.0 x 10 1 


2.0 x 10" 


-3 


990123 


14.4 


5.3x10" 


-4 


5.1 x 10 4 


2.1 x 10" 


-3 


990510 


11.9 


3.1x10" 


-4 


3.2 x 10 1 


1.2x10" 


2 


991208 


8.8 


1.0x10" 


-2 


2.0xl0~ 2 


1.0 x 10" 


-1 


991216 


8.5 


1.7 x 10 


-4 


3.0X10- 2 


7.0 x 10" 


-2 


000301c 


10.8 


2.0x10" 


-2 


7.0xl0~ 3 


1.2x10- 


■ 1 


000418 


10.9 


5.0x10- 


-3 


7.9X10- 1 


4.5 x 10" 


-2 


000926 


9.5 


3.5x10- 


-2 


1.5X10" 1 


1.4 x 10" 


-1 



(25) 



s = 
s = 2 



(26) 



In Table 2, we provide the theoretical estimations of the 
magnitude of the reverse shock emission for the ten bursts 
in our sample (for a homogeneous external medium and as- 
suming the RS parameters to be same as the FS) , to deter- 
mine if these bursts would have had a bright optical flash. 
These results were obtained numerically using an accurate 
calculation of the cooling and the self-absorption frequen- 
cies, which can also be found in Table 2, and the flux is 
found to be consistent with the analytical estimate given 
above. Our fluxes are somewhat smaller than reported in 
Soderberg & Ramirez- Ruiz (2002). The difference is perhaps 
because the RS falls in a regime that is neither Newtonian 
nor relativistic where the usual asymptotic approximations 
are not very accurate, and one needs a more accurate calcu- 
lation for this intermediate case (Nakar & Piran 2004 have 
made a similar point). 

We see that for six of the ten bursts the optical flux 
is between 8.5 and 11-th magnitude, and the cooling fre- 
quency at deceleration is small (less than or of order the 
R-band frequency). With the cooling frequency either be- 
low or dropping below the optical rather quickly, and no 
electrons left in the RS capable of radiating in the optical 
band, we would not see the expected ~ l/t 2 falloff, but a 



a Reverse shock R-band magnitude at deceleration with pa- 
rameters determined from afterglow modeling, 
k Redshift not known for this burst, z = 1 used. 



more rapid falloff of ~ l/i A (Kumar & Panaitescu, 2000). 
This falloff is fast, but even so, some of these bright optical 
transients could be seen for a few hundred seconds by rapid 
followup observations with a limiting magnitude of R ~ 15. 

The remaining four have much fainter optical emission 
and large cooling frequency. The magnetic field parameter 
determined from afterglow modeling for these bursts is much 
lower than those with bright reverse shock emission. With 
the high cooling frequency (at least a factor of 10 above the 
observing band frequency of 2 eV), we expect this emission 
to exhibit t~ 2 falloff; however, this faint emission may be 
hidden by brighter forward shock afterglow emission. 

For those bursts which are fit equally well in an s — 
2 medium (970508, 000418, 991208, 000301c, 991216), the 
reverse shock flux at deceleration has been calculated using 
the afterglow parameters in PK02. We find that the flux is 
typically a few times larger in the s=2 model than uniform 
ISM case discussed above. 

According to table 2, the peak R magnitude for the 
reverse shock emission for 990123 is 14.4 (using the forward 
shock values for all parameters in the reverse shock) whereas 
the observed peak flux was R = 8.9 (Akerlof et al. 1999). 
This suggests, as has been pointed out by Zhang et al. 2003, 

that eb was larger in the RS by a factor of about 10 2 than in 

3/2 

the FS. Since v c oc e B , by making ebt ~ 0.07, the cooling 
frequency has been lowered from 51 keV to 100 eV. The R 
band is at about 2 eV, so the cooling frequency is still well 
above the optical at deceleration, allowing for the ~ l/t 2 
falloff that was observed. 

Bursts with small optical flux at deceleration have small 
eb and/or e' e in the RS and their cooling frequency is gener- 
ally high, allowing the optical emission falling off as ~ l/t 2 
to occur. If eb in the RS for these bursts were larger than 
the FS, as found for 990123 (Zhang et al. 2003) & 021211 
(KP03) , the flux will be boosted to levels comparable to that 
of the other brighter bursts. 

Besides 990123 and 021211 there have been no obser- 
vations of a bright and quickly fading early afterglow for 
any other bursts. There have been some early (time since 
onset of burst < 0.01 day or ~ 10 minutes) detections in 
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the optical, for example 030329 and 021004, but the emis- 
sion was not falling off as ~ l/t 2 . There have been about 
18 bursts with upper limits published in the GCN Circulars 
(Barthelmy et al. 1995) and in the literature (e.g. Akerlof et 
al. 2000, Kehoe et al. 2001), for emission between the GRB 
time and 0.01 day (9 of these bursts reported in the GCNs 
had later optical afterglow detections). Searching the GCN 
Circulars with the GRBlog website (Quimby et al. 2003), 
we find that the burst upper limits range from R ~ 10 for 
030115 (Castro- Tirado et al. GCN 1826) at early times to 
R ~ 20 for XRF 030723 (Smith et al. GCN 2338) closer 
to 0.01 day. None of the bursts from the sample in this pa- 
per have upper limits available. However, if the bursts with 
available upper limits are representative of the total GRB 
population, then it is possible that the bursts in our sam- 
ple would have had similar upper limits, i.e., roughly 14-15th 
magnitude at ~ 500s. So there is a disagreement between the 
theoretical expectation and the observational upper limit. 

There are several possible resolutions for this apparent 
discrepancy. The small optical flux could be due to much 
smaller magnetic field in the RS compared with the FS; for 
990123 and 021211, however, eb in the RS was inferred to 
be larger than the FS, which perhaps might not be the com- 
mon situation. Another possibility is that the deceleration 
time of GRB fireball is of order an hour instead of the burst 
duration of a few tens of seconds (assumed for the GRBs 
considered in this paper). Since the peak optical flux from 
RS is proportional to ~ t^ 1 (see eauations l25l fc l2T>I . the flux 
will be reduced by ~3-5 magnitudes and therefore consistent 
with the observational upper limits of Kehoe et al. (2001). 
In this case the early lightcurve should be rising, and subse- 
quently turn over to a steep decay at the deceleration time. 
The most likely explanation for a typically faint optical flux 
in our view is related to the low cooling frequency in the 
RS. We see from the table 2 that cases with a bright opti- 
cal flash have cooling frequency below the optical band at 
deceleration, in which case the lightcurve should decline as 
~ t~ 3 and fade below the detection limit of 14-15 magnitude 
in a few hundred seconds. 

The observational situation (whether bright reverse 
shock emission is typical or not) should become clearer when 
Swift is launched in September 2004. 



5 DISCUSSION 

We have explored the possibility that prompt 7-ray emission 
for a selected sample of 10 long duration GRBs might arise 
in the external shocks. These bursts had good multiwave- 
length afterglow data and temporal coverage which enabled 
Panaitescu and Kumar (2002) to determine their energy, 
jet opening angle, density of the surrounding medium, and 
microphysics parameters for the shock. We compared the 
observed peak flux and the peak frequency of the time av- 
eraged spectrum during the 7-ray burst with the theoretical 
extrapolation of the afterglow emission to the middle of the 
burst duration. 

The motivation for considering external shocks for the 
generation of 7-ray emission for these 10 bursts is that most 
of these bursts are not highly variable, which is the primary 
reason for invoking internal shocks. Moreover, the efficiency 
for the production of 7-rays for these bursts is found to be 



very high - in excess of 50% - which is difficult to under- 
stand in the internal shock models. In the particular case of 
GRB 990123 the 7-ray pulse width was of the same order as 
the time when the prompt optical emission from the reverse 
shock peaked i.e. the deceleration time. This means that the 
radius where 7-rays were generated was roughly the same as 
the deceleration radius, thereby suggesting a forward shock 
origin for gamma-rays. Moreover, the low energy spectral 
slope a (/„ oc v a ) for 990123 was 0.4, an observation that 
is difficult to understand in the internal shock models which 
have generally very low cooling frequency and therefore have 
a = -0.5. 

We find that it is not possible that the forward 
shock synchrotron afterglow emission, extrapolated back to 
prompt GRB duration, can explain the flux and peak fre- 
quency of nine of these ten bursts; only in the case of 970508, 
which was a single pulse FRED burst, does the extrapolation 
of the afterglow match up with the 7-ray emission property. 
Moreover, it turns out that even when we take esf (the 
energy fraction in magnetic field) and ISM density during 
the gamma-ray burst to be completely arbitrary, instead of 
the same as determined from afterglow modeling, we still 
cannot reconcile the gamma-ray observations with the theo- 
retical calculations for seven of the ten bursts in our sample 
in the forward shock synchrotron emission model. The rea- 
son for this is that the forward shock is highly radiative at 
early times i.e., the cooling frequency during the gamma-ray 
burst is smaller than the synchrotron injection frequency, 
and therefore the flux at the peak of the vf v spectrum is 
independent of the density of the surrounding medium and 
EBf, which happen to be the parameters with large uncer- 
tainty in afterglow modeling. 

Two of the bursts in the sample, 980519 & 000418, can 
be understood as synchrotron emission in the forward shock 
provided that esf ~ 1 during the burst , a value that is larger 
by a factor~ 10 2 than what we find during afterglow at 
~ lday; it is unclear if this is a physically sensible solution. 

We have also considered inverse-Compton scattering of 
synchrotron emission from reverse or forward shock off of 
material in the forward and reverse shock regions, and find 
that it is not possible to explain the 7-ray emission, except 
possibly for 980703, with a reasonable set of parameters. In 
particular, the only solutions we found are when the circum- 
burst medium is taken to have a wind like density profile 
with the density parameter about hundred times larger than 
the density of a typical Wolf-Rayet star wind and a few 
order of magnitude larger than the density determined from 
afterglow modeling. 

Adding large density clumps to the external medium 
only increases the peak flux and frequency of vf v by a sig- 
nificant amount when the density of the clump and ejecta 
are similar when they collide to produce the 7-ray emis- 
sion, in which case the early afterglow emission is extremely 
bright and hard to miss. Also, adding e pairs to the ejecta 
decreases the inverse Compton flux, unless the electron dis- 
tribution power law index is <2, for any synchrotron-IC scat- 
tering scenario considered. Thus, neither of these two possi- 
bilities seem likely to explain the 7-ray emission properties 
of the bursts in our sample. 

So, we are therefore forced to conclude that there must 
be another way to explain the GRB emission. The widely 
accepted internal shocks model might be the solution. How- 
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ever, considering the problem of efficiency, special cases of 
FREDs, and the problems for the internal shocks for 990123 
discussed in §1, we feel it is prudent to explore other possi- 
ble mechanisms such as the conversion of Poynting flux to 
radiation. 

As for the reverse shock emission, we find that at least 
50% of these ten bursts have bright prompt optical flashes, 
9- 11th mag, provided that the shock parameters - the en- 
ergy fraction in electrons and magnetic field - in the reverse 
shock are the same or larger than the value in the forward 
shock. However, five of these bursts have cooling frequencies 
below the optical band and therefore the RS flux will decline 
very steeply with time (roughly as t~ 3 ) and could easily go 
undetected after a few deceleration times. Those bursts with 
dimmer early emission generally have high cooling frequen- 
cies, and are assumed to exhibit the expected l/t 2 fall off. 
Although they are dim, they may still be detected. 

If this sample of 10 bursts is representative of long du- 
ration GRBs then we expect very bright, rapidly fading 
(~ i -3 ), prompt optical flashes accompanying many 7-ray 
bursts. The RS emission is particularly bright just above the 
synchrotron self-absorption frequency of ~ 10 _1 eV where 
we expect the observed flux to be of order a Jansky, and 
declining rapidly with time since the cooling frequency is 
typically of the same order as the absorption frequency. It is 
also possible that the deceleration time is much longer than 
the GRB duration, which would reduce the predicted optical 
flux from the reverse shock. In this case we would expect to 
observe a dim, rising early optical afterglow lightcurve, turn- 
ing over to a steep descent (t~ 2 or t~ 3 ) after the deceleration 
time. 

These issues will be resolved in the Swift era when we 
will have excellent early time coverage in the optical band 
for a few hours for many bursts. Future measurements of 
the early afterglow should enable us to determine if bright 
optical emission from the reverse shock is common, and thus 
determine the nature of the explosion, i.e. whether the ex- 
plosion is baryonic, leptonic, or Poynting-flux dominated. 
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